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ABSTRACT 

We present the results of recent Chandra, XMM-Newton, and Hubble Space Telescope observations of the 
radio-loud (RL), broad absorption line (BAL) quasar PG 1004-1-130. We compare our new observations to 
archival X-ray and UV data, creating the most comprehensive, high signal-to-noise, multi-epoch, spectral 
monitoring campaign of a RL BAL quasar to date. We probe for variability of the X-ray absorption, the 
UV BAL, and the X-ray jet, on month-year timescales. The X-ray absorber has a low column density of 
Ah = 8 X 10^*^-4 X 10^' cm“^ when it is assumed to be fully covering the X-ray emitting region, and its 
properties do not vary significantly between the 4 observations. This suggests the observed absorption is not 
related to the typical “shielding gas” commonly invoked in BAL quasar models, but is likely due to material 
further from the central black hole. In contrast, the C IV BAL shows strong variability. The equivalent 
width (EW) in 2014 is EW = 11.24±0.56 A, showing a fractional increase of /SEW/{EW) = 1.16±0.11 
from the 2003 observation, 3183 days earlier in the rest-frame. This places PG 1004H-130 among the most 
highly variable BAL quasars. By combining Chandra observations we create an exposure 2.5 x deeper 
than studied previously, with which to investigate the nature of the X-ray jet and extended diffuse X-ray 
emission. An X-ray knot, likely with a synchrotron origin, is detected in the radio jet 8" (30 kpc) from the 
central X-ray source with a spatial extent of 4" (15 kpc). No similar X-ray counterpart to the counterjet 
is detected. Asymmetric, non-thermal diffuse X-ray emission, likely due to inverse Compton scattering of 
Cosmic Microwave Background photons, is also detected. 
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1. INTRODUCTION 

Active Galactic Nuclei (AGN) can significantly affect the 
properties of their host galaxies through mechanical feedback 
provided by winds and jets. Past (v ^ 0.01 -0.1c) outflow¬ 
ing winds manifest themselves observationally as blueshifted, 
broad absorption lines ( BALs) in the UV/optic al spectra of 
^ 15% of quasars (e.g.. iHewett & Poltzl l2003h . although it 
is likely that they exist in most quasars and we only de¬ 
tect them when our line-of-si ght passes directly th rough 
the outflowing material (e.g., IWevmann et al.l 119911) . As 
this material is not static, variability is expected in the 
BALs on a range of timescales characteristic of, e.g., vari¬ 
ations in the ionization parameter and/or changes in the 
flow structure or disk rotation. There have been many stud¬ 
ies o f the UV variability of BALs in s amples of quasars 


(e.g.. 

lBarlowlll993l: iLundgren et al.ll2007 

iGibson et al.l|2008 

201C 

; iCanelluno et al 1120 11 L 201 2l 120 1 31: 

Filiz Aketal.112012 

2013 

; IWildv et all 12014 IWelling et all 12014. These And 


that BAL variability is both more common and stronger on 
longer rest-frame timescales and in weaker and higher veloc¬ 
ity BALs, with the dom inant parameter bein g the equivalent 
width (EW) of the BAL (iFiliz Ak etalJlT^ . 

In the disk-wind scenario, the wind is launched from 
the accretion disk at ~ lOOORg and is driven by UV line 
pressure. For a central source emitting typical levels of 
X-rays, this requires some “shielding gas” to be located 
between it and the wind to prevent over-ionization of the 


outflow (e.g., [Murray et al.l 119951 : iProga et ^ l2000h . This 
gas may be responsible for the absorption commonly ob¬ 
served in the X-ray spectra of BAL quasars which is well 
fit with either a partially covering and/or a partially ionized 
absorption model with typical column densities of Ah ^ 
10^1-” cm-^ (e . g..lGallagheretalJ[2()()l iGiustini et alJUOOl 

iFan et al.]l2009HStreblvanska et alfcoiC ). Thi s X-ray absorp¬ 

tion is also observed to vary in some case s (iGallagher et al.l 
l2004t IChustini et aDl201 IbISaez et al.l[2012l) . possibly as a re¬ 
sult of complex dynamics related to rotating, infalling, or 
outflowing material. Compared with non-BAL quasars of 
a similar UV lu minosity, BAL qu a sars are underluminous 
in X-rays (e.g ., iGrem et al.l 119951 iGreen & Mathuii 119961: 

IGallagher et alJlT999ll2006tlBrandt et alJl200C)l) . but, after cor¬ 

recting for the effects of absorption, a majority are found 
to have similar underlying X-ray spectral propert i es to those 
of no n-BAL quasars (e.g., IGallagher et al.l l2002t iPage et akl 
120051) . However, s ome BAL quasar s are likely to be intrin¬ 
sically X-ray weak (iTeng et alJ 2014t iLuo et al1l2014l) and in 
such cases, shielding gas is not required to ensure launching 
of the wind. 

Approximately 10% of AGN have powerful, highly col¬ 
limated and relativistic jets which emerge from the in¬ 
ner regions (^ 10 Rg) of the AGN and emit strongly 
at radio and X-ray frequenci es via synchrotron and Syn¬ 
chrotron Self Compton (SSC; [Band & Grindlavl 119861) pro¬ 
cesses. Radio-loud quasars (RLQs) which possess such 
jets have higher X-ray luminosities than radio-quiet quasars 
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(RQQs; e.g.. IZamorani et aL 119811 : IWorrall et alJ Il987t 
iMiller et alJ 120111 : IScottetal, 120111) . an d flatter (harder) 
X-ray power-law spec tral indices (e.g., iReeves & Turner! 
[20001: LSaez et all (2?^ as a result of this additional iet - 
linked X-ray emission component (iSambruna et al.l Il999h . 
X-ray counterparts to exten ded radio jets have b een ob¬ 
serve d in many AGN (e.g ., ISambruna et al.l 12004 IWorralll 
l2009t iMarshall et al.l 120111) . Diffuse X-ray emission has 
also been detected tra cing the direction of r adio jets and un¬ 
seen counteriets (e . g.. ISchwartz et al.ll20()^ and radio lobes 
(e.g.. lCroston et al]l2005l) and is likely due to inverse Comp¬ 
ton scattering of Cosmic Microwave Background (CMB) pho¬ 
tons (IC/CMB) by relativistic jet electrons. 

It was originally thought that q uasars could not ho st both 
a BAL wind and a radio jet (e.g., IStocke et alJll992h . How¬ 
ever, many subsequent discoveries of radio-loud (RL ) BAL 
quasars have now been made (e.g.,[Becker et al.ll 19971120001: 

[Brotherton et al.ll998l:lMenou et alJ200ll) although BALs are 

still not found in quasars with t he highest levels of radio¬ 
loudness (e.g.. iGregg et al.ll2006t IShankar et ^l2008l) . This 
suggests that the radio jet and BAL wind may affect each 
other, despite being emitted from different locations within 
the AGN. It is perhaps more likely that physical properties 
of the central engine such as mass accretion rate and the 
spin of the black hole are responsible for the presence of ei¬ 
ther phenomenon, although no significant differences in the 
black hole masses or the Eddington ratios of radio-quiet (RQ ) 
and RL BAL quasars have been found (iBruni et alJ 120141) . 
RL an d RQ BAL quasars also show similar o ptical prop- 
erties (iRunnoe et'ffl 120131; iRochais et alJ l2014h . However 
iBaskin et akl (1201 ^ suggest that differences in the BAL prop¬ 
erties of RQ and RL quasars are due to a softer UV/optical 
SED, indicated by a lower He II EW, which are preferen¬ 
tially found in RL BALs. BAL varia bility in RLQs doe s 
not depend strongly on radio loudness (lEiliz Ak et alJl20f3h . 
but lobe-dominated RLQs may show gr eater fractional BAL 
variability than core-dominated RLQs dwelling et al.ll20l4) . 
When considering variability on the same timescales, RLQs 
show a lower fractio nal variability than RQQs, ^ 40 ± 20% 
dwelling et al.ll201^ . RL BAL quasars are X-ray weak com¬ 
pared to RL non-BAL quasars as expected, but to a lesser de¬ 
gree than the difference between RQ BAL and RQ non-BAL 
quasars. This suggests that the X-ray a bsorption does not af¬ 
fect all of the jet-linked X-ray emission dMiller et al.ll2009l) . 

PG 1004H-130 (also known as PKS 1004H-13 and 4C 13.41) 
is one of the best-studied exa mples of a RL BAL quasa r. 
It is optically bright {B = 15.8: iVer on-Cettv & Veronl 20101) . 
low redshift (z = 0.2406, lEracleous & HalpernI 12004 ). and 
has a black hole mass of Mbh = 1-87;!;q45 x 10‘^Mq 
dVestergaard & PetersonI |2006|) . International Ultraviolet 
Explorer (lUE) and Hubble Space Telescope (HST) spec¬ 
tra of PG 1004H-130 show blueshifted absorption (v ^ 
10,000 kr ns~Q in the high- ionization lines O VI, N V, Si IV, 
and C IV d Wills et al.ll 19991). and it has a C IV BALnicity in- 
dejQ of BI ^ 850 kms~^ dWills et al^ll999^ . It also shows 
stron g C IV BAL variabi lity, with large fractional changes in 
EW dwelling et akllMl) . 

Like most BAL quasars, PG 1004H-130 is X-ray weak 


/-3000 


/25000 

0.9. 


where /(v) is the normalized flux as a function of velocity displacement from 
the rest-frame wavelength of the emission line. The BI does not include the 
first 2000 kms“^ of the absorption trough iWevmann et al.lll99in . 


and was undetected i n an early Einstein observation 
dElvis & Eabbianol 11984) . The first X-ray detection and 
spectra of PG 1004H-130 were obtained with XMM-Newton 
in 2003 and Cha ndra in 2005, and were reported in 
IMiller et"^ d2006l: hereafter M06). Although the XMM- 
Newton spectrum did not show any significant X-ray absorp¬ 
tion, the Chandra spectrum obtained 20-months (observed- 
frame) later, revealed complex X-ray absorption best-fit with 
a partial-covering model (Ah = 1 . 2 ;|;og 4 xlO^^ cm“^, / = 
0A9^j,l, r = L37!°j^). PG 1004H-130 has also been ob- 
serve d by the Nuclear Spe ctroscopic Telescope ARray (NuS- 
TAR; iHarrison et al.l 120131) . It was weakly detected in the 
10-20 keV band, suggesting that the shielding X-ray gas 
is either Compton thick or that PG 1004H-130 is intrinsically 
X-ray weak, alt hough neither sc enario was strongly preferred 
over the other (iLuo et al.ll2013l) . If the absorbing material is 
Compton thick, we would expect to detect a strong Ee Kg 
emission line, but this is not observed (M06 : ILuo et al.ll2013l) . 
Its absence could be due to X-ray em ission related to a jet di¬ 
luting the spectrum (iLuo et al .||20 13|), a scenario also favo red 
by broadband modeling (iKunert-B^raszewska et anil2015l) . 

PG 100 4-H130 is RL with a radio loudness parameter of 
R* = 209 (I Wills et all 119991: where R* = fscm/fpop\, and 
R* > 10 indicates a RLQ) and its exte nded radio emission is of 
a HY MOR (HYbrid MORphology; iGonal-Krishna & Wiital 
I 2 OOOI) nature. The surface-brightness of the radio lobe to the 
SE dec reases outwards i.e., a E anaroff-Riley I (ERI) classifi¬ 
cation (lEanaroff & Rilevll 19741) . while the surface-brightness 
of the NW lobe increases outwards (ERII). This asymmetry is 
likely due to twin jets propagating into dissimilar large-scale 
environments, i.e., the density of the surrounding medium is 
higher to the SE so that the je t is more quickly decollimated 
(iGonal-Krishna & Wiital l2000l) . The X-ray imaging analysis 
reported in M06 revealed an X-ray counterpart to the radio jet 
extending to the SE ^ 8" from the nucleus and upstream from 
the peak of the radio jet emission. Large-scale diffuse X-ray 
emission was also detected, extending 40-50" outward from 
the nucleus with a flux of ^ (4.2-4.5) x 10“^^ ergs“'cm“^ 
(assuming a power-law model with L = 1.8; M06). 

In this work, we present results from newly obtained X-ray 
and UV observations of PG 1004H-130, and compare them 
to archival data. This allows us to probe for variability of 
the X-ray absorption, the C IV broad absorption line, and the 
X-ray jet emission on month-to-year timescales. By com¬ 
bining multiple Chandra observations we obtain the deepest 
imaging to date of an X-ray jet in a HYMOR radio source, 
with 2.5 X greater sensitivity than the original data used in 
M06. This allows us to place further constraints on the spatial 
and spectral properties of the X-ray jet, and larger-scale dif¬ 
fuse X-ray emission. This study of PG 1004H-130 represents 
the most comprehensive, high signal-to-noise (S/N), multi- 
epoch, X-ray and UV spectral monitoring campaign of a RL 
BAL quasar, with results applicable to this class of quasars 
as a whole. A description of our new and archival X-ray and 
UV observations, and details of the data reduction are given in 
Section|2] In Section|3]we present an analysis of the X-ray and 
UV emission from the nuclear source, including the X-ray ab¬ 
sorption properties and the C IV BAL variability. In Section|4] 
we present analysis of the extra-nuclear X-ray jet and diffuse 
emission. We discuss our results in Section |5| and summarize 
our conclusions in Section|6] A standard cosmology in which 
Ho = 69.7 kms“'Mpc“*, Ua = 0.7185, and Hm = 0.2815 is 
assumed throughout (iHinshaw et al.ll20T3l) . This gives a lumi¬ 
nosity distance to PG 1004H-130 of dv= 1216 Mpc, and hence 











































































































Scott et al. 


3 


an angular scale of 3.8 kpc per arcsec. Errors quoted on X-ray 
quantities are 90%, unless otherwise stated. 

2. OBSERVATIONS 

2.1. Chandra Data Reduction 

PG 1004H-130 has been observed t wice with the Advan ced 
CCD Imaging S pectrometer (ACIS: iGarmire et al.l120031) on 
board Chandra dWeisskopf et al.ll2000t) . We present results 
from both a recently obtained 60 ks observation (taken in 
2014 Mar) and a 40 ks observation taken in 2005 Jan which 
was previously reported in M06. Details of both observations 
are given in Table|T] We reduce each d ata set in the same, con - 
sistent way using CIAO version 4.6.1 dFruscione et al.ll2006ll 
and version 4.6.2 of the CALDB. The chandra_repro 
script is used to apply the latest calibration data. Neither 
observation suffers from significant background flaring, but 
we use a standard 3tT clipping of outliers from the 2.4-6 keV 
light curve to select good-time intervals (GTIs) where the 
count rate was < 0.23 s“' or < 0.21 s“' in 2005 and 2014, 
respectively. Improved screening for potential cosmic-ray 
background events was carried out on the 2014 data which 
were taken in VFAINT mode. Neither observation suffers 
from pile-up in the nuclear source. We use specextract 
to obtain spectra for the point-like nucleus using a circular 
extraction region with a 4" radius centered on the X-ray 
source position, and background spectra are extracted from 
a 16" radius circle, offset from the source position. These 
spectra are binned to a minimum of 20 counts per bin to 
ensure the validity of statistics during the subsequent 
spectral fitting. Auxiliary response hies are generated 
using the latest calibrations to include the time-dependent 
quantum-efficiency (QE) degradation due to contamination 
build-up on the ACIS optical-blocking hlter. X-ray spectra of 
the jet emission are also extracted from a rectangular region 
of size 6.4" x 4.3", offset by 8" from the nuclear source, 
and at a position angl^l of 130° from the source position 
in both observations (see Section 0. Background spectra 
are extracted from the same 16" radius circular regions 
as used previously, but in this case no binning is applied, 
an d the spectra are modeled in XSPEC using the Cash statis¬ 
tic (ICashll979h due to the low number of counts in the spectra. 


2.2. XMM-Newton Data Reduction 

We also consider two XMM-Newton (iJansen et al.l 2001) 
European Photon Im aging Camera (EPIC: lTurner et akl 2001 ; 
IStriider et al.l 1200 Ih observations of PG 1004H-130; a more 
recent 60 ks observation taken in 2013 Nov as part of AO-12 
and a 20 ks archival observation from 2003 May, hrst re¬ 
ported in M06. Neither of these observations were affected by 
substantial background Haring, therefore GTIs were dehned 
as times where the count rate was < 0.4s“' in the pn detector 
and < 0.35 s“' in the MOS detectors (although for the 2013 
observation this was reduced to 0.115s“' to exclude 2 short 
Hares). Final exposure times and details of the observations 
are given in Table [T] Data from each observation were 
reduced in a consistent and standard way using SAS version 
13.5.0. Spectra of the point-like nucleus were extracted 
from a circular region with a radius of 24.6", as in M06, 
centered on the source position. Neither observation suffers 
from pile-up of the nuclear source. Circular regions on the 

^ All angles are measured from North, through East. 


same CCD but offset from the source with radii of 43" and 
34" were used to extract background spectra from the MOS 
and pn detectors, respectively. Auxiliary response Hies and 
redistribution matrices were generated using arfgen and 
rmf gen. The resulting spectra from the two MOS detectors 
were combined, and both the resulting MOS and pn spectra 
are binned to a minimum of 20 counts per bin using the Ftool 
grpphaQ 


2.3. HSTData Reduction 

We used th e Co smic Origins Spectrograph 
(COS; iGreen et^ 1201 2^ onboard the HST to obtain 
a current-epoch Near-UV (NUV) spectrum covering 
Aobs = 1670-2120 A which includes the C IV BAL region 
of PG 1004H-130. A 2562 s exposure was obtained with 
the G230L grating at a central wavelength of 2950 A on 
2014 Jan 20. The G230L grating was selected to provide 
spectral resolution exceeding archival UV observations while 
spanning the entire C IV emission line and BAL absorption 
region on the A stripe. The target was centered in the 
cross-dispersion and along-dispersion directions using an 
ACQ/PEAKXD and ACQ/PEAKD acquisition sequence. 
Four focal plane grating offset positions (FPPOS=l^) were 
employed to accumulate the spectrum on different regions 
of the detector. Data were taken in the TIME-TAG mode to 
retain the arrival time information for each event. Standard 
pipelin43 processing steps were applied, including Doppler 
correction for the shift induced by the orbital motion of HST, 
non-linearity correction for detector deadtimes at high count 
rates, Hat-Held correction for pixel sensitivities, wavelength 
calibration with a flashed line-lamp, bad pixel flagging, 
extraction of the background-subtracted one-dimensional 
spectrum, conversion from counts to flux, and conversion 
to heliocentric wavelengths. The final spectrum has a S/N 
ratio of ^ 10 across the C IV BAL. No signiflcant temporal 
or spectral variability is present over the duration of the 
observation. 


3. X-RAY AND UV ABSORPTION VARIABILITY 
3.1. Long-term Light Curve 

Figure [T] shows long-term light curves for PG 1004H-130. 
The top panel shows the observed 0.5 - 8 keV X-ray fluxes 
(not corrected for Galactic or intrinsic Ah), estimated from 
an absorbed power-law model in which the power-law slope 
is free to vary. XMM-Newton fluxes are shown in red, and 
Chandra fluxes are shown in blue, both with 90% errors. 
Shown in black is the flux extrapolated from a NuSTAR ob¬ 
servation in 2012 Oct (iLuo et al.ll2013h . the green arrow in- 
dicates a 3(7 upper limit from an early Einstein observation 
(lElvis & Fabbianoll 1984ft . and the magenta arrow is an upper 
limit derived from the ROSAT All Sky Survey (IVoges et al.l 
119991) . Each was extrapolated to a 0.5-8 keV flux by assum¬ 
ing a power-law slope of F = 1.5, typical of radio-loud AGN 
(e.g.. lReeves & Turneill200dl) . 

The bottom panel shows historical monitoring giving an 
indication of the typical optical variability (^ 0.5 mag) of 
PG 1004H-130. V-band magnit udes are est i mated from the 
unflltered P magnitudes given in lSmith et al.1 (119931) . differen- 

^ littp://heasarc.gsfc.nasa.gov/docs/software/ftools/ 

^ littp://www.stsci.edu/hst/cos/documents/handbooks/datahandbook/ 
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Table 1 

X-ray Observation Log 


Start Time 

Obsid 

Observatory 

Detector(s) 

Exposure Time 

Nuclear Source Counts'^ 

Observed 0.5-8 keV Elux* 

(UTD) 




(ks) 


(x 10“l3ergs“lcm^3) 

2003 May 04 19:06:59 

0140550601 

XMM-Newton 

EPIC-MOS 

21.51 

1171)35 

40+0.17 




EPIC-pn 

18.08 

1574)j‘ 


2005 Jan 05 16:44:27 

5606 

Chandra 

ACIS-S 

39.29 

1763)j3 

4 S7+0.29 
^•o/_0.30 

2013 Nov 05 05:58:10 

0728980201 

XMM-Newton 

EPIC-MOS 

63.14 

2399)50 

7 -70+0. 10 

^•'°-o.io 




EPIC-pn 

54.99 

3598)01 


2014 Mai- 20 10:42:47 

16034 

Chandra 

ACIS-S 

60.97 

1599)^1 

Q4-b0.22 

■^•■^^-0.21 


Notes. 

‘^The number of observed counts in the nuclear source extraction region. Vox XMM-Newton this includes counts in the energy range 0.5 — 10 keV, from 
a circular extraction region with a radius of 24.6^^; for Chandra this includes counts in the range 0.5 - 8 keV from a A" radius region. 

^’The obseiwed 0.5 - 8 keV fluxes are determined from an intrinsically absorbed power-law model. For the XMM-Newton obseiwations the MOS and pn 
spectra are modeled simultaneously giving a single flux estimate. 


tial V-band photometry given in iGarcia et alJ ( r999h. differ- 
ential /?-band photo metry given in IStalin et alJ (l2004ll . pho¬ 
tometry from SDSS (lYork et ^l2000h . and unfiltered CCD 
magnitudes from the Catalina Real-Time Transient Survey 
(iDrake et al.ll200^ . 

There exists a well-known correlation be tween the X-ray- 
to-op tical power-law slope parameter, q;o}< 0 (iTananbaum et al.l 
Il979h . and the rest-frame UV luminosity, due to the physical 
intei~p lay between the UV disk photons and the X-ray corona 
(e.g. lyignali et alJ l2003t IStrateya et alJ 120051 : ISteffen efaP 
l20()6t lYoung et alJ 120101: iLusso et al] 1201 Oh . Adopting 
log / 2500 A“ 30.78 from Ishen et al.l (1201 Ih . and after cor- 
recting for the ex cess X-ray luminosity expected in RLQ 
(iMiller et al.l 1201 ih . an X-ray flux of Tb. 5 - 8 key = (75-95) x 
10“^^ ergs“^cm“^ is expected. This is 15-35x higher than we 
obserye in any of the 4 XMM-Newton and Chandra obserya- 
tions showing the persistent X-ray weakness of PG 1004H-130. 

Flux estimates of the same object measured by the EPIC 
and ACIS detectors are known to show discrepancies with 
ACIS measuring fluxes that a r e typically "^10- 15% higher 
(e.g., iNeyalainen et alJ 1201 Ot iTsuiimoto et al.l 1201 Ih . The 
2005 ACIS flux is ~ 40% higher than the 2003 flux measured 
by EPIC, and the 2014 ACIS flux is ^ 20% higher than 
the 2013 EPIC flux. These differences, which are larger 
than those typically expected, are suggestiye of a real flux 
change in the source between the obseryations, rather than 
simply calibration differences. When we consider the two 
XMM-Newton, or the two Chandra obseryations separately, 
which reduces calibration effects, we see a significant flux 
decrease on timescales of ~ 8 years in the rest-frame of 
PG 1004H-130. A similarly significant decrease is obseryed in 
the hard band flux (2-8 keV). The lower flux measured by 
NuSTAR is also suggestiye of a real flux change in the source. 


3.2. X-ray Spectral Fitting of the Nuclear Source 

The grouped X-ray spectra were fit with physically moti- 
yated models using XSPEC yersion 12.8 (lArnaudI 19961) . Each 
spectral model includes absorption fixed at the yalue of the 
Galactic component in the directio n of the source . Thi s 
is estimated using equation 7 from IW illingale et al.| (|2 013h . 
Ahi = 3.70 X 10^*’ cm“^ from iDickey & LockmanI (I1990I) . and 

^ «ox = 0.381og (/x//b) where/x and/o are tlie rest-frame flux densities 
at 2keV and 2500 A, respectively. 

^ 4.5(7 for XMM-Newton and 6.1(7 for Chandra. 
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Figure 1. Long-term light curves of PG 1004-1-130. Top - the observed 
0.5 - 8 keV fluxes estimated from an absorbed power-law model. Fluxes 
from both XMM-Newton obseiwations considered in this work are shown in 
red, and the fluxes from both Chandra observations are shown in blue. Each 
is plotted with a 90% error bar. The observed 0.5 - 8 keV flux extra polated 
from a NuSTAR obseiwation is shown in black <Luo et af] 120131) . The 
green arrow represents a 3cr upper limit from an early Einstein obseiwation 
(Elvis & Fabbiandll984l) . and the m agenta arrow repre sents an upper limit 
from the ROSAT All Sky Survey iVoges et^ 1 19991). B ottom - Filled 
orange sy mbols show optical m agnitudes fromISmit h et alJ ( 19931 : upwards 
triangles), IGarcia e t alj 119991: squ ares), IStalin et*^ <20041: downwards 
triangles), SPSS fVork et all 12006: circle), and the Catalina Real-Time 
Transient Survey <Drake et alJ2Q09l : diamonds). 


E(B — V) = 0.0331 mag from ISchlafly & EinkbeineJ (1201 Ih . 
resulting in N^^’ = 4.1 x lO^** cm~^. Absorption i s mod eled 
using tbabs with abundances fromlWilms et aP (120001) and 
cross-sections from IVerner & Eerlandl (119961) . Unless stated 
otherwise, the XMM-Newton spectra were modeled oyer the 
energy range 0.5 - 10 keV. The MOS and pn data were fit si¬ 
multaneously using the same parameters with a freely yarying 
constant added to the mod el to account for cal ibration offsets 
between the two detectors (iMateos et alJl200^ . The Chandra 
ACIS spectra were modeled oyer the energy range 0.5-8 keV. 

Each of the 4 data sets was first modeled with a simple 
power-law fit to the high-energy data (aboye 2 keV), absorbed 
only by a Galactic contribution. Eigure|2] shows the unfolded 
spectra (top) and the ratio of the data to this power-law model 
(bottom). This clearly shows a flat power-law slope, as ex¬ 
pected for radio-loud AGN such as PG 1004-(-130. When the 
hard-band power-law fit is extrapolated down to 0.5 keV, the 
presence of intrinsic absorption at lower energies is also clear, 
and appears to be larger in the Chandra spectra than in the 
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XMM-Newto n spectra. Such absor ption is expected for BAL 
quasars (e.g. lGallagher et ^l2002h . 



Figure 2. X-ray spectra from each of the 4 epochs. The XMM-Newton 
MOS and pn spectra are fit simultaneously (data from MOSl and MOS2 
are combined) and are plotted in the same color on this figure for clarity. 
Each of the data sets are modeled with a power law (with only a contribution 
from Galactic absorption) over 2 — 8 keV {Chandra) and 2-10 keV 
{XMM-Newton). This model is shown extrapolated down to 0.5 keV. The 
vertical grey line indicates the observed energy at which the neutral iron 
emission line at a rest-frame energy of 6.4 keV is expected. This feature is 
visible in the 2013 XMM-Newton spectra. The top panel shows the unfolded 
spectra with E x Ef{E) plotted on the y-axis such that a power-law slope 
of r = 2 would be horizontal in the figure. The bottom panels shows the 
ratio of the data to the power-law model, with the residuals at lower energies 
clearly indicating the absorption present in both Chandra spectra, and the 
2013 XMM-Newton spectra. 


We therefore fit each of the data sets with an absorbed 
power-law, XSPEC model tbabs*ztbabs*po, i.e., the ab¬ 
sorber at the redshift of the source is both fully neutral and 
fully covering the source. The best-fitting parameters for 
each model are listed in Table |2] Confidence contours (at 
1, 2, & 3cr) from this absorbed power-law fit are shown in 
Figure [3l where different colors and linestyles are used for 
each of the 4 data sets. Crosses mark the location of the 
fit parameters for which is a minimum, although none 
is shown for the 2003 XMM-Newton data as Ah is uncon¬ 
strained in this fit. The results from the two Chandra ob¬ 
servations (2005, shown in dashed blue, and 2014, shown 
in solid green) are highly consistent, showing constrained ab¬ 
sorption column densities of Ah, 2005 = (O.lOfg Qg) x 10^^ cm“^ 
and Ah, 2014 = ^ cm“^, respectively. However, 

there are clear differences between these Chandra data and 
the results from both XMM-Newton observations, which show 
steeper photon indices and lower levels of absorption. 

The X-ray detectors of ACIS and EPIC are known to yield 
different spectral fitting results even during simultaneous ob¬ 
servations of the same object, or i n observations of non - 
variable objects (e.g., galaxy clusters: [Nevalainen et al.ll2010l) . 
By comparison with the typical F and Ah differences found 
in simultaneou s ACIS and EPIC observations of the BE Lac 
PKS 2155-304 (llshida et al.ll2011h . we conclude that the 2013 
XMM-Newton results are likely to be consistent with both the 



0.8 1.0 1.2 1.4 1.6 1.8 

Power-low index, P 


Figures. Confidence contours (1, 2 & 3ij), for an absorbed power-law fit 
(tbabs*ztbabs*po) for each of the 4 data-sets. Crosses indicate the 
fit parameters which give the minimum value of The Chandra data 
appears to be consistent between each epoch. The 2003 XMM-Newton data 
yields significantly different spectral parameters; differences that appear to 
be larger than those expected from known calibration differences between 
the detectors. The 2013 XMM-Newton data also shows differences, but of a 
smaller size, such that they may be attributed to calibration effects. 


2014 and 2005 Chandra observations after accounting for 
these cross-calibration differences. Only the original 2003 
XMM-Newton observation appears to show PG 10043-130 
with different spectral parameters (F, Ah). This was previ¬ 
ously observed and discussed in M06. 

However, the detectability of absorption is dependent on 
the spectral quality. As the 2003 XMM-Newton spectra con¬ 
tain fewer counts than are required to significantly detect ab¬ 
sorption of Ah 0.08 X10^^ cm“^ (the level meas ured in the 
2013 XMM-Newton data), 100% of the time (see IScott et al.l 
l2012l for detectability curves), this may be limiting our ability 
to significantly detect absorption even though it is present in 
the source. If instead, we fix F to the value obtained from 
the fit to high energies (> 2 keV), as shown in Figure |2l we 
do measure a significant absorption component with Ah 2003 = 
(0.08!^^) X10® cm This is consistent with the absorption 
measured in the 2013 XMM-Newton data using the same spec¬ 
tral modeling. We therefore suggest that the spectral proper¬ 
ties of PG 1 004 h- 130 indicated in both XMM-Newton observa¬ 
tions may be consistent, and therefore the amount of intrinsic 
absorption has not varied dramatically between all 4 observa¬ 
tions. 

If each of the data sets are modeled with an absorbed power- 
law in which F is fixed to the values obtained from the fit 
to high energies (> 2 keV), steeper power-law slopes and 
higher absorption columns are obtained (see Table |2|i. How¬ 
ever, we not e that the F values rema in slightly flatter than ex¬ 
pected (e.g.. lReeves & Turneill200^ . which may indicate that 
the shape of the underlying continuum above 2 keV is still 
affected by the absorption. If the spectra are modeled with 
a power-law using only energies greater than 4 keV, slightly 
steeper spectral indices are obtained but they remain formally 
consistent (< Icr) with the previous estimates. Adopting these 
steeper F values in the subsequent absorption modeling leads 
to a possible increase in the absorption column density of a 
factor < 2. 

Thus far the spectral modeling has only considered absorp¬ 
tion due to neutral material that is fully covering the cen¬ 
tral source. While this is useful for a simple and consistent 
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Table 2 

X-ray Spectral Fitting Results 


Data 

Model 

r 

Ah 

/ 

log? 

HO 

F-test 




(X10'^'^ cwr^) 


[erg s~^ cm] 

(%) 

(%) 

XMM-Newton 

Power law 

1 T^+O.OS 
-0.05 




17.3 


2003 

PL + absorption 

1 3A+0.07 
i-jO_oo5 

<0.06 



15.7 

0 


PL (2-10 keV) 

-0.14 




82.2 



PL + absorption 

1.50 (fix) 

0.08^:“ 



7.7 



PL + partially covering absorber 

1.50 (fix) 

4.96fi'i 

A 23+0.07 
^•^•^-0.08 


18.4 

99.7 


PL + ionized absorption 

1.50 (fix) 

0.58^:“ 


2-OK.n 

12.7 

97.6 

Chandra 

Power law 

^•^^-0.05 




8.3 


2005 

PL + absorption 

^•^^-0.10 

A iq-fO.IO 
^•^^-0.08 



36.5 

99.9 


PL (2-8 keV) 

1 42+0.20 
^ -^^-0.20 




19.8 



PL + absorption 

1.42 (fix) 

0 25+0.07 
^•”^^-0.06 



12.8 



PL + partially covering absorber 

1.42 (fix) 

1 40+0.84 
^ -^^-0.59 

^•°^-0.07 


36.9 

99.8 


PL + ionized absorption 

1.42 (fix) 

0.55tfi 


0.89^]^ 

18.7 

94.4 

XMM-Newton 

Power law 

1 iQ-fO.03 
^•^^-0.03 




27.6 


2013 

PL + absorption 

1 97-f0.05 
^•^^-0.05 

0.08^:“ 



38.7 

99.5 


PL (2-10 keV) 

1 'J':3+0.09 

i-JJ-0.09 




44.5 



PL + absorption 

1.33 (fix) 

^•^^-0.03 



34.2 



PL + partially covering absorber 

1.33 (fix) 

0 92+^-28 
^•^^-0.66 

'^•■^^-0.08 


40.0 

96.3 


PL + ionized absorption 

1.33 (fix) 

^•^°-0.21 


1 42+0.66 
-^^-2.94 

36.4 

87.2 

Chandra 

Power law 

^ -^^-0.07 




45.1 


2014 

PL + absorption 

^•^^-0.12 

0 17+0.13 

U.i /-O 12 



63.3 

98.9 


PL (2-8 keV) 

1 '^S+0.21 




93.1 



PL + absorption 

1.38 (fix) 

^••^°-0.09 



36.4 



PL + partially covering absorber 

1.38 (fix) 

^••^'-0.77 

A 54+0.14 
-0.09 


62.2 

99.7 


PL + ionized absorption 

1.38 (fix) 

0 71+0.88 

^•'^-0.30 


1 97+0.55 
^•^^-0.64 

51.2 

97.8 


Notes. The spectral fits are carried out using data in the energy range 0.5 — 10 keV for XMM-Newton and 0.5 — 8 keV for 
Chandra, unless stated otherwise. Parameters are quoted with 90% errors, g = L/nr^ with units of ergs^' cm; log g= —3 for 
neuhal material, and log +6 for fully ionized material. HO values greater than 1% indicate the model is an acceptable fit 
to the data. F-test values greater than 99% indicate that the model including the additional parameter is statistically preferred. 
Despite differences in the fitting methodology, the best-fitting parameters for the partially covering absorber model for the 
2005 Chandra data are consistent with those given by M06. 


quantification of the amount of absorption present in each 
observation, the presence of more complex X-ray absorption 
is th ought to be common in both RQ and RL BAL qu asars 
(e.g.. lGallagher et ^120021 l2006t iBrotherton et al.ll2005h . We 
therefore also fit the spectra with a partially covering absorp¬ 
tion model (tbabs* zpcf abs*po in XSPEC) and a partially 
ionized absorption model (t babs*zxipcf*po), in which 
ZXIPCF (iReeves et al.ll2008h is implemented with / fixed at 
1.0 to model a partially ionized, but fully covering, absorber. 
We fix the F value to that found for high energies to ensure 
constrained Na, ^ and / parameters. These are listed in Ta¬ 
ble 121 along with null-hypothesis probabilities (HO), which 
if greater than 1% indicate the model is an acceptable fit to 
the data, and F-test probabilities, for which we say values 
greater than 99% indicate that the model including the addi¬ 
tional component is statistically preferred. 

In all 4 observations, the F-test probability and HO values 
are higher for the partially covering absorption model than the 


partially ionized absorption model, although both provide ac¬ 
ceptable fits to the data. For both models the best-fitting spec¬ 
tral parameters are consistent when comparing the 2 Chan¬ 
dra observations or the 2 XMM-Newton observations. This 
again suggests that the properties of the X-ray absorption have 
not changed significantly on long timescales 8 yrs in the 
source rest-frame). 

We test for the presence of iron emission lines in each of 
our data sets by adding Gaussians to a power-law fit at high 
energies, i.e., the model tbabs* (po+zgauss) . We use 
the Cash statistic and unbinned spectra as the increased bin 
size in the grouped data could limit the detection of a narrow 
line. We fix the width of the line to u = 0.01 keV and fix the 
rest-frame energy of the line to the values 6.4 (for neutral 
iron), 6.7 (He -like), and 6.9 k eV (H-like). In agreement 
with M06 and iLuo et all (1201 3h . we find no strong evidence 
for iron emission lines at these energies. We obtain only 
upper limits on the EW in all cases except for the 6.4 keV 
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line in the 2013 XMM-Newton spectra. This is unsurprising 
as these spectra have the highest spectral resolution and 
highest number of counts of all those we consider. The EW 
measured is EW 6.4 kev = eV, and the line can be seen 

in the residuals shown in Eigure|3 We also consider a joint 
fitting of all 6 X-ray spectra (MOS, pn, and ACIS) and obtain 
EWe 4 kev = 90±40 eV (consistent with the upper limit of 
EW 6.4 keV < 105 eV obtained by M06), EWe.? keV < 70 eV, 
and EWe.g keV < 70 eV, suggesting that a low level of neutral 
iron emission is present. 


3.3. Broad Absorption Line Variability 

In this section we compare our newly obtained NUV spec¬ 
trum from HST COS to archival spectra taken by lUE in 1982 
and 19 86, and a HST Space Telescope Imaging Spectrograph 
(STISi lWoodgate et alJl998h spectrum taken in 2003, in order 
to look for variability of the C IV BAL. Details of all 4 obser¬ 
vations are given in Table [3 We also include the S/N of the 
spectra and the observed UV continuum fluxes, measured be¬ 
tween Aobs = 1750-1825 A (Arest =1411-1471 A). With the 
exception of a ^ 2 (t flux decrease between the 2003 and 2014 
HST spectra, there does not appear to be significant variation 
in the UV flux. 

Each of the 4 spectra were modeled with a Voigt profile 
and a linear function to represent the C IV emission line and 
the local continuum, respectively. These models are shown 
by the dashed lines in the top panel of EigureH] The spectra 
were then divided by the best-fitting model, resulting in the 
ratio spectra shown in the bottom panel of EigureH) 

This shows a mini-BAL with a EWHM of 1553 kms“^ 
centered on Arest 1542 A. The offset of the absorp¬ 
tion from the expected emission wavelength of C IV (A = 
1551 A) corresponds to a typical outflow velocity of 
Tout ~ 1745 kms“' for the mini-BAL. The BAL has a EWHM 
of 3372 kms“^ (in the 2014 spectrum), indicating a large 
range of velocities observed along the line-of-sight within the 
accelerated outflow. These range from Vout 3700 kms“' to 
Tout 11800 kms“' ('^ 0.04c), with the velocity of the deep¬ 
est part of the BAL trough being at Vout ^ 8000 kms“\ 

The EW of the BAL and mini-BAL are directly measured 
from the ratio spectra. The edges of the BAL region are de¬ 
fined as the wavelengths at which the BAL reaches a level 
which is 90 % that of the continuum following the standard 
definition of IWevmann et al.l (1199 Ih . In the 2014 spectrum 
this corresponds to a range of Arest = 1491 - 1532 A for the 
BAL, and Arest = 1537-1550 A for the mini-BAL. The EW is 
determined simply as a sum of the flux in all pixels within this 
region, and the error estimate includes contributions from the 
continuum placement error, statistical spectral noise, and the 
uncertainty in the region used. The EW of the BAL and mini- 
BAL in the newly obtained 2014 spectrum are 11.24 ± 0.56 A, 
and 3.38 ±0.42 A, respectively. These are listed in Table [3 
along with EW estimates for the 3 archival spectra determined 
by the same methocfl. The final rows of Table [3 compare the 
EW values of the BAL measured in each of the spectra, with 
the fractional change in EW being given by 

AEW ^ EWnew-EWold 

{EW) 0.5(£Wnew + £Wold)' ^ ’ 

^ We use a differ ent normalization and modeling of the spectra to that of 
IWelling et al.l resulting in slightly different EW values. 


Eigure |4] (bottom) shows the C IV BALs for each of the 
4 spectra, smoothed with a Gaussian of 8 A width. This 
clearly shows the strongly variable C IV absorption that 
PG 1004-1-130 is known to possess, i.e., the BAL strength in¬ 
creases from the 1982 lUE observation (blue) to the 1986 lUE 
observation (red), before diminishing greatly in depth in the 
2003 HST observation (green). Eurther, dramatic variability 
is revealed by the latest 2014 HST observation (black), which 
shows the BAL has deepened. In addition to large changes 
in EW, the shape of the BAL trough also varies significantly 
between the observations. In comparison, the mini-BAL at 
a smaller blueshift shows little variability in both strength or 
shape, with the exception of the 2003 spectrum in which the 
trough is much shallower. 

The 2014 HST spectrum does not cover the wavelength of 
other possible BAL troughs (e.g., O VI, N V), nor do we see 
evidence for a Si IV BAL. 



1450 1500 1550 1600 

Rest-frome wavelength (A) 


Figure 4. UV spectra of PG 1004+130, covering the C IV 
emission/absorption-line region. Archival lUE spectra from 1982 and 
1986 are shown in blue and red respectively, and an archival HST STIS 
spectrum from 2003 is shown in green. Our recent HST COS spectrum (2014 
Jan) is shown in black. The top panel shows each of the spectra normalized 
to the same scale. The dashed lines show the best-fitting model for each 
spectrum, composed of a linear component to model the local continuum, 
and a Voigt profile for the C IV emission line. The bottom panel shows the 
resulting ratio spectra after each spectrum has been divided by its best-fitting 
model. They have been smoothed with a Gaussian of 8 A width. The regions 
over which the EW of the BAL and mini-BAL are determined are shown by 
the vertical, black, dotted lines. 


4. EXTRA-NUCLEAR X-RAY EMISSION 

The 2005 Chandra observation of PG 1004-1-130 revealed 
both X-ray emission ~ 8 " SW of the nucleus, thought to be 
related to the ERI radio jet, and diffuse X-ray emission sur¬ 
rounding the source and extending out to ~ 40-50" (M06). 
In this work, we combine the data from both Chandra obser¬ 
vations giving an increased exposure time of ~ 100 ks (2.5x 
deeper than previously) with which to further investigate these 
features. 
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Table 3 

UV Observations and BAL Parameters 



Observation 1 

Observation 2 

Observation 3 

Observation 4 

Date 

1982-04-24 

1986-01-12 

2003-03-31 

2014-01-20 

MJD 

45083 

46442 

52729 

56678 

Mission 

lUE 

WE 

HST 

HST 

Instrument 



STIS 

COS 

Grating 



G230L 

G230L 

Central wavelength (A) 



2376 

2950 

Resolution (A) 

6 

6 

~5 

1 

Exposure time (s) 

22260 

18900 

2292 

2562 

UV flux (X10“*^ erg s“' cm“^) 

9.29 ±0.87 

8.03 ± 1.01 

8.53±0.41 

7.20 ±0.49 

Signal-to-noise (S/N) 

10.6 

7.7 

20.6 

14.6 

Mini-BAL EW (A) 

3.38 ±0.32 

4.14±0.38 

2.74 ±0.25 

3.38 ±0.42 

BAL EW (A) 

3.03 ±0.56 

13.38 ±0.64 

2.98 ±0.34 

11.24±0.56 

Comparison of BAL to 2014: 

T (days) 

9346 

8251 

3183 


AEW (A) 

8.21 ±0.79 

-2.14±0.85 

8.26 ±0.66 


(EW) (A) 

7.14 ±0.40 

12.31 ±0.43 

7.11 ±0.34 


AEW/{EW) 

1.15±0.13 

-0.17 ±0.07 

1.16±0.11 


Comparison of BAL to 2003: 

T (days) 

6163 

5068 



AEW (A) 

-0.05 ±0.66 

-10.40 ±0.72 



(EW) (A) 

3.01 ±0.33 

8.18±0.36 



AEW/(EW) 

-0.02 ±0.26 

-1.27±0.10 



Comparison of BAL to 1986: 

r (days) 

1095 




AEW (A) 

10.35 ±0.85 




(EW) (A) 

8.21 ±0.43 




AEW/{EW) 

1.26±0.12 





Notes. The UV flux listed is the observed flux between Aobs = 1750 - 1825 A 
(Arest = 1411 — 1471 A) located bluewards of the C IV BAL. The bottom rows compai'e 
the BAL properties measured in each of the spectra. r is the time between observa¬ 
tions in the rest-frame of PG 1004+130 and is calculated as r = (MJDnew — MJDoid)/(l + 
z). AEW = EWnew - EWold is the change in the EW, and (EW) = 0.5(£Wnew + EWom) 
is the average EW. The ratio of these two quantities gives the fractional change listed. 


4.1. X-ray Jet Emission 

Figure|5](top left) shows the 0.5-4 keVQX-ray image pro¬ 
duced from the addition of the 2005 and 2014 Chandra obser¬ 
vations. The 2014 data were reprojected onto the co-ordinate 
system of the 2005 data to account for a 0.7" offset between 
the source positions. The spectral extraction regions used for 
the nuclear source and the X-ray jet emission are indicated by 
the green circle and rectangle, respectively. The blu e contours 
are pr oduced from 4.9 GHz (6 cm) VLA radio data (iFomalonll 
119811) with a resolution of 1.5", and show flux density levels 
of 0.75 and 1.25 mJy beam“*. Figure |5] (bottom left) shows 
the same X-ray image, but on a larger scale. In thi s case, radio 
data from the FIRST survey (iBecker et al.lll995h at 1.4 GHz 
(21 cm) with a resolution of 5.4" are used to give contours at 
1.5, 3, 8, and 15 mJy beam“* (overlaid in blue). This clearly 


shows the HYMOR (iGopaTKrishna & Wiitall2000l) nature of 
the radio emission, with the SE jet showing a decreasing 
surface-brightness (FRI), whereas the NW counterjet shows 
an increase in surface brightness at larger radii from the nu¬ 
cleus, typical of the FRII class. X-ray emission is observed 
following the same direction as the SE radio jet, albeit at 
closer radial distances than the peak of the radio emission. 

Eigure |5] (top right) shows the central region of the com¬ 
bined 0.5-4 keV X-ray image which has been rebinned to 
a subpixel size of 0.05" and adaptively smoothed using a 
Gaussian kernel to achieve a minimum S/N of 3. X-ray 
contours corresponding to 0.0005, 0.001, 0.002, 0.005, 0.01, 
0.02, and 0.05 counts per square pixel are overlaid in green. 
These correspond to flux densities between 6.6 x 10“^^ and 
6.6 X 10“^^ergs“'cm“^Hz“^ arcsec“^, assuming a photon en¬ 
ergy of 1 keV, for which the ACIS QE is ~ 0.5 and the ef- 


There are few counts at energies greater than 4 keV. 
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fective area is ^ 400 cm^0 This image again clearly shows 
the X-ray knot which lies upstream of the radio jet emission 
at a distance of ^ 8 " (projected distance of 30 kpc) and posi¬ 
tion angle of 130° from the nuclear source. No similar X-ray 
emission is observed in relation to the NW counterjet. 

Spectra of the X-ray jet emission to the SE of the nuclear 
source are extracted from each individual Chandra observa¬ 
tion using the rectangular region shown on Figure |5] (top 
left). There are 29 X-ray counts (0.5-8 keV) in both the 
2005 and 2014 jet spectra, ~ 2-3 of which are likely back¬ 
ground counts, with so me d ue to the larger-scale diffuse X-ray 
emission (see Section I4.21 l. The spectra are modeled with 
a power law absorbed only by a Galactic contribution. The 
2005 spectrum is relatively flat with T = 1. 63^0 57 , whereas 
the 2014 spectrum is comparatively steep, T = 2.21^^^^. How¬ 
ever, due to the large (90%) errors, this variation is not sig¬ 
nificant (0.9cr). We also simultaneously model both spec¬ 
tra to obtain a better constraint on the power-law slope; T = 
1.96;!^ 52 - ^ decrease in the flux of the jet is measured, 
falling from Eo.5-8kev,2005 = (5.47+^;^^) x 10“'^ergs“'cm“2 to 
Tb, 5 - 8 kev, 20 i 4 = (2.86;|;[ [ 2 ) X 10“'^ergs“'cm“^. This variation 
is largely due to the difference in the best-htting model as 
when r is fixed to 1.96 for each observation, the measured 
fluxes differ by only ^ Icr and show a decrease of ^ 30%. 

Radial surface-brightness prohles for the SE X-ray jet emis¬ 
sion were created from the 0.5-4 keV X-ray images which 
were hrst rebinned to a sub-pixel size of 0.05". Concentric 
annular sectors covering a 30° arc from 115- 145°, with a 
thickness of 1" were used to determine the jet prohle. Fig- 
ure| 6 ](top left) compares the jet prohles from the 2005 (black) 
and 2014 (red) observations. This shows that the distribution 
of the X-ray emission has not varied signihcantly between the 
2 epochs. An increase in the X-ray surface brightness is clear 
at ^ 8 " from the nuclear source, corresponding to the peak 
of the X-ray jet emission. This knot is also clearly extended 
by 4". Figure | 6 ] (top right) shows the jet prohle from the 
combined 2005 and 2014 data (black). This is compared to a 
prohle created from the remaining 330° of the annular sectors 
(red). This clearly shows the excess X-ray surface brightness 
to the SE, which is signihcantly higher than the background 
level (estimated from a 8 " x 8 " box and indicated by the 
dashed blue line). 


4.2. Diffuse X-ray Emission 

Figure |5] (bottom right) shows the wider-scale X-ray im¬ 
age, smoothed in the manner described in the previous sec¬ 
tion. This image shows more clearly the low-level diffuse 
X-ray emission surrounding PG 1004H-130 which was previ¬ 
ously reported by M06. The emission largely appears to fol¬ 
low the direction of the SE FRI jet and the NW FRJI coun¬ 
terjet, rather than forming a spherical halo around the nuclear 
source. The diffuse X-ray emission to the SE appears to trace 
the full extent of the lobe, and extends out to ^ 65", corre¬ 
sponding to a projected distance of ^ 250 kpc. The emis¬ 
sion to the NW extends out a similar distance as to the SE 
(^ 55" 210 kpc), but does not trace the full extent of the 

FRII radio lobe. The X-ray emission to the NW is also weaker 
than to the SE. There are 129;!;22, 0.5-4 keV background- 
subtracted counts in a 50° wedge to the SE, extending from 

^http://cxc.harvard.edu/proposer/POG/html/chap6.html 
(Chandra handbook, chapter 6) 


16" to 63" (marked on Figure |5] bottom right, in greenf*^ 
and 86 ;!; 2 i counts in a similar wedge along the NW counter 
jet. The contributions from point sources have been excluded 
and the values are corrected for any subsequent loss of area. 
Owing to the increased exposure time of the combined obser¬ 
vations, and the subsequent increase in counts, we are able to 
conduct a basic spectral analysis of the diffuse emission. We 
extract spectra from the SE and NW wedge regions and bin 
them to a minimum of 15 counts per bin using grppha. A 
power-law model gives E = 1.65 ^q7?, jv = 19.6/13 and a 
flux of Eo.5-4keV = (1.02/[’j^) X 10“®erg S^'cm”^ (Eo.5-4keV = 
1-81/o 44 X 10"^^ergs”') for the SE emission and E = 1.95/[ oI, 
= 9.9/10 and fb, 5 - 4 kev = (4.82/j g^) x 10”'^ergs”'cm“^ 
(L(), 5 - 4 kev = 8 . 56/3 2 ^ 10"** ergs”*) for the NW emission. This 
conhrms the lower flux of the NW emission, and gives consis¬ 
tent power-law slopes. Although we cannot formally reject a 
bremsstrahlung model (it gives similar reduced values and 
HO > 1%), the kT value obtained is unconstrained, tending to 
higher temperatures and thus a power-law distribution. 

Figure |6] (bottom) shows the radial surface-brightness 
prohles of the diffuse X-ray emission. The numbers of counts 
were extracted from 15 concentric arcs with equal thickness 
between a radius of 0-60". The prohle of the SE diffuse 
emission (shown in black) was extracted from a 50° arc 
between angles 95° and 145°, the diffuse emission along the 
NW counterjet (shown in red) was extracted from a 50° arc 
between angles 275° and 325°, and remaining portions of 
the arcs are used for the comparison prohle (shown in green). 
The background level (shown by the blue, dashed line) was 
estimated from a 30" x 30" box at a distance > 80" from 
the central source position. This hgure shows that the X-ray 
surface brightness along the jet and counterjet directions is 
higher than that in the remaining directions. The surface 
brightness of the SE diffuse emission is also generally higher 
than that of the NW emission, although not significantly at 
each radius. In general, both profiles show a smooth decrease 
with increasing radial distance from the source. 

5. DISCUSSION 

5.1. The Nature of the X-ray and UV Absorption 

The X-ray absorption commonly observed in BAL quasars 
may be due to “shielding gas” located close to the central 
black hol e which prevents the over-ionization of the BAL 
wind ( e.g..lMurrav et alJl995tlProga et alJ200CilGibson et al.l 
l2009t IWu et al.r 20101) . Quasar-wind simulations suggest 
this gas should be varia b le on timescales of months-to- 
yrs (e.g., ISim et all 120101 l2012h . X-ray absorption vari¬ 
ability has been observed in some BAL quasars such as 
PG 2112-1-059 whic h showed a dramatic i ncrease in column 
density over 3 yrs (iGallagher et al.l l2004ll . and PG 1126- 
041, a mini-BAL quasar which showed rapid variability on 
timescales of months -to-hours due to highly ionized absorbers 
(iGiustini et al.ll20fl]) . However, i n the sample of 11 BAL 
quasars studied systematically by ISaez et alT (120121) . excep¬ 
tional X-ray variability (e.g., uncovering events where the 
shielding gas largely moves out of the line-of-sight), com¬ 
pared to non-BAL quasars, was not seen. 

In this work we have shown that the X-ray absorption prop¬ 
erties of PG 1004-1-130 measured in 4 observations taken up 

Note that the SE region used here differs slightly from that used by 
M06, as the X-ray emission appears to be diametrically opposed to the NW 
emission in this deeper exposure. 
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Figure 5. Top left - 0.5-4 keV X-ray image from the combined 2005 and 2014 Chandra observations. The X-ray emitting nuclear source and the X-ray knot 
in the SE jet are clearly visible, and the extraction regions used for the nuclear source (circle) and the jet (rectangle) are shown in green. Radio contours from 
the VLA at 4.9 GHz (6 cm) corresponding to 0.75 and 1.25 mJy bearn^* are overlaid in blue. Top right — An adaptively smoothed version of the top left image. 
Shown in green are X-ray contours corresponding to count rates of 0.0005, 0.001, 0.002, 0.005, 0.01, 0.02, and 0.05 counts per square pixel. Blue contours show 
the same VLA data as in the top left image. Bottom left - The same 0.5 — 4 keV X-ray image shown in the top left image, but including the extended regions 
of the system. Blue contours show lower angular resolution radio data from FIRST (at 1.4 GHz/21 cm) at 1.5, 3, 8, and 15 mJy bearn^*. Bottom right - An 
adaptively smoothed version of the bottom left image. The black dots to the SE and W are unrelated X-ray sources, likely background AGN. The green wedges 
con'espond to those used to determine the amount of diffuse X-ray emission along the paths of the SE jet and NW counterjet. The blue contours show the same 
FIRST data as in the bottom left image. At the redshift of PG 1004-1-130, 15” corresponds to a projected physical distance of 57 kpc. 


to 8.8 yrs (rest-frame) apart do not differ substantially, regard¬ 
less of the spectral model used. When a fully covering and 
neutral absorber is assumed, the column densities range from 
A^h = 8 X 10^°-4 X 10^^ c m~^, which are among the lowest val¬ 
ues typically observed ((Ga llagher et al.| l200a iGiustini et al.l 
l2008t[^n et alJ2()0^IStreblvanska et al.l20ldl) . However, the 
power-law slopes obtained in the modeling are F < 1.5, per¬ 
haps indicating that the true levels of absorption in th e source 
may be higher by a factor of ~ 2 (see Section 13.21) . In the 
case of a partially covering absorber, the column densities ob¬ 
tained are naturally higher, A^h 10 ^^ cm“^, and more consis¬ 
tent with those for other BAL quasars. The X-ray absorption 
can be well modeled with either a partially covering, or an 
ionized, absorber (see Table|2]) and the spectral parameters of 
these models (Ah, f, and 0 are also consistent for each obser¬ 
vation. 

As the X-ray absorption properties do not appear to vary, 
the absorption is perhaps more likely to be due to material 
located further from the central source, rather than the typi¬ 


cal shielding gas expected on small scales. In the case that 
PG 1004-1-130 is intrinsically X-ray weak, such a shield is 
not required to ensure launching of the BAL wind. Alter¬ 
natively, if PG 1004-1-130 is intrinsically X-ray normal, the 
shielding gas is required to be highly Compton-thick (Ah ~ 
7 X 10^"* cm~^: iLuo et £0120 1 3h . Absorption with column den¬ 
sities this high would not be detectable in Chandra or XMM- 
Newton spectra as all X-rays below ^10 keV would be ab¬ 
sorbed. In both scenarios the X-ray emission we observe is 
therefore predominantly from the jet; either because the con¬ 
tinuum source is weak in comparison, or because we only see 
X-rays from the portion of the jet extending beyond the shield¬ 
ing gas. This X-ray emission may then be absorbed by mate¬ 
rial further from the black hole, naturally explaining the lack 
of variability and lower column densities observed. However, 
we cannot rule out variability occurring at a level below that 
which we are able to signihcantly detect with our spectra. A 
column density increase of more than a factor of 2 would be 
detected as a 3(7 change between the Chandra spectra. 
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Figure 6. Radial surface brightness profiles. In each case the profiles are extracted from the 0.5 — 4 ke V Chandra ima ges which were rebinned to a sub-pixel 
bin size of 0.05". Surface brightnesses are calculated as the number of counts per square pixel and the lGehr^ 119861) formulation of errors is used. Top left 
- A comparison of the SE jet profiles for the 2005 data (black) and the 2014 data (red), extracted from concentric 30° arcs at angles of 115— 145°, centered 
on the source position, and with a thickness of 1". The blue dashed line gives an indication of the background level (~ 0.02 counts per square pixel in 2005, 
~ 0.03 counts per square pixel in 2014) which is determined from an 8" X 8" box. Top right — The SE jet profile extracted from the combined 2005 and 2014 
Chandra observations (black), using the same regions as in the top left panel. The profile shown in red is extracted from the remaining 330° annular sectors for 
comparison. The blue dashed line indicates the background level (~ 0.05 counts per square pixel in the combined data). Bottom - Profiles of the diffuse X-ray 
emission extracted from concentric 50° arcs to the SE (black) and NW (red). These are compared to a profile extracted from the remaining 260° from each 
concentric ring (green). Fifteen equal-thickness rings are used from 0 - 60" radial distance from the central source position. The blue, dashed line indicates the 
background level (0.02 counts per square pixel) determined from a 30" X 30" box. 


In the scenario in which PG 1004-1-130 is intrinsically 
X-ray normal and is absorbed by Compton-thick shielding 
gas, a strong Fe Ka emission line at 6.4 keV with an EW 
of 1000-2000 eV is expected (e.g.. iGhisellini et'H] 119941: 
iMatt et al]ll996h . Although we do measure a low level of 
neutral iron emission in a joint fittin g of all 4 X-ray spectra 
(EWg qjjev = 90±40eV; see Section 1X2] ). it is not as strong 
as might be expected. This may be du e to the dilution of the 
spectrum by jet-linked X-ray emiss ion (iLuo et al.ll2013h . 

The X-ray light curve in Section [3T| shows significant flux 
variability between the 4 X-ray observations. Given the con¬ 
stant level of X-ray absorption present, these changes are 
more likely related to intrinsic variations within the contin¬ 
uum source itself (in the case that PG 1004H-130 is intrinsi¬ 
cally X-ray weak), or fluctuations in the jet. 

In stark contrast to the X-ray absorption properties, the 
C IV BAL in PG 1004H-130 shows large-amplitude variabil¬ 
ity. The archival observations we compare with our newly 
obtained HST COS spectrum probe rest-frame timescales 
of « 3-26 yrs during which large absolute and fractional 


changes in EW are observed (see Table O. However, the 
largest changes are between the 2003 and 1986 spectra, 
i.e., A£'W = -10.4±0.72Aand A£'W/(£'W) =-1.27±0.10. 
The magnitudes of these changes place PG 100 4-H130 among 
the m ost variable RE and RQ BAL quasars (e.g.. lPiliz Ak et all 
120131) . Our measurements probe long rest-frame timescales, 
on which BAL variability is observed to be stronger (see Sec¬ 
tion [T]). In Eigure^we compare the absolute and fractional 
BAL variability of PG 1004H-130 (r ed, solid points) to the sam¬ 
ple of RL BAL quasars presented in IWelling et al.l (l2014t open 
black circles). This highlights the large variability of 
PG 1004H-130 in relation to other RL BAL quasars, even de¬ 
spite the longer rest-frame timescales probed by our observa¬ 
tions. 

Variability in BAL troughs is likely caused by 2 main ef¬ 
fects; clumps of material in the outflow moving across our 
line-of-sight to the continuum source, or the material respond¬ 
ing to changes in the incident ionizing radiation. When only 
portions of the BAL trough spanning a narrow velocity range 
show changes, perhaps with some portions strengthening 
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Figure?. The C IV BAL variability in PG 1004+130, and other RL BAL 
quasars. Filled red squares show the values for PG 1004+130 and open 
black circles s how v alues from the sample of RL BAL quasars presented in 
IWelling et all j2014i) — excluding PG 1004+130. The top panel shows the 
relationship between the absolute variability and rest-frame time interval. 
The bottom panel shows the fractional variability. 

while others weaken (e.g. JGibson et al.ll2008tlCapellupo et all 
l2012h . clouds at different velocities, and hence distances from 
the continuum source, entering and leaving the line-of-sight 
is a plausible explanation, although i t could also be a resul t 
of velocity-dependent optical depths (ICapellupo et al.ll2012h . 
However, many observations of BALs show co-ordinated 
changes in separate troughs of the same ion, for example, 
distinc t C IV BAL troughs tend to strengthen or weaken to¬ 
gether (iFiliz Ak et al.l20l'3h . In this case, variations occurring 
over such a large velocity range are unlikely to be the result 
of co-ordinated movements within the outflow wh ich would 
disperse the structure (e.g.. iRogerson et"^ 1201 Ih . Instead, 
these, and global changes in which the enti re profile of a sin - 
gle wide BAL trough varies uniformly (e.g.. lGrier et al.l2015h . 
are likely a result of a change in the ionizing EUV continuum. 
However this can only occur if the absorption trough is not 
too highly saturated. Portions of the trough with low optical 
depths may show larger variations as they are more suscepti¬ 
ble to changes in ionization, than those with only moderate 
optical depths. 

The changes observed in the BAL profile of PG 1004H-130 
suggest a combination of these scenarios. This is most ap¬ 
parent when considering the BAL trough in the 1986 observa¬ 
tion, which appears to be comprised of a main trough centered 
on Arest ~ 1511 A (vout ^ 7830 kms“*) and an additional ab¬ 
sorption feature with a smaller blueshift superimposed. This 
component, at Arest 1526 A (vout ~ 4870 kms“'), does not 
appear in the other 3 observations (there is a hint of a sim¬ 
ilar component in the 1982 spectrum but its depth is within 
10 % of the continuum level and could therefore be noise), and 
its narrow width suggests it may be due to cloud of material 
in the line-of-sight during the observation. Making the as¬ 
sumptions that such a cloud is spherical, moving perpendicu¬ 
lar to the line-of-sight and passing in front of a negligible size 
source, we can place an upper limit on the radius of the cloud 
of r < 1.3 X 10*^ cm < 0 .04 pc, as it must have traveled at 
least its diameter, at v ^ 4870 kms“', in the 6163 day period. 

The remainder of the BAL trough in the 1986 spectrum 
has a similar shape to that in 1982, although both the width 
and depth have increased by more than a factor of 2. Such 
a monolithic change is more likely to be the result of the 
entire outflow varying in response to a change in the ion¬ 
izing continuum. However, the UV flux does not change 


significantly between these observations. Similarly, there is 
no significant flux difference between the 1986 and 2003 
observations, despite a large decrease in the BAL strength, 
AEW = -10.40±0.72 A. The flux decrease between the 2003 
and 2014 observations is only significant at 2cr but does corre¬ 
spond to a strengthening of the BAL (AEW = 8.26 ± 0.66 A). 
The correlation between the observed UV flux and BAL EW 
may not be as strong as expected as the flux is measured 
over Arest = 1411- 1471 A. While this should give an indi¬ 
cation of the UV flux variability, it is not directly measuring 
the unobservable EUV range which would be responsible for 
any changes within the outflow. Previous observations have 
shown the EUV fl ux can vary more than the UV flux at longer 
wavelengths (e.g.. lMarshall et aT]ll997^ . 

If a change in the ionizing continuum is the main con¬ 
tributor to the BAL variability, we might expect to see co¬ 
ordinated variations in the mini-BAL. This feature is centered 
on Arest ~ 1542 A corresponding to an outflow velocity of 
Tout ~ 1745 kms“' which is over 10000 kms“' slower than 
parts of the BAL trough. Eigure |4] shows that the mini-BAL 
in the 2003 HST spectrum is significantly weaker than in the 
other 3 spectra. This corresponds well to the BAL which is 
also extremely weak in this spectrum. However, in general, 
this feature shows much less variability than the BAL; in par¬ 
ticular there is no corresponding deepening in 1986 or 2014 
when the BAL is considerably stronger. This is likely due 
to the feature being more highly saturated, and therefore less 
susceptible to changes in the ionizing continuum, but may 
also suggest that material is indeed physically distinct from 
the higher velocity BAL and likely located much at a much 
greater distance from the central source. 

The 2014 BAL trough shows a similar shape to the main 
BAL trough in the 1986 spectrum, but is offset in wavelength. 
This could be due to an acceleration of the entire outflow 
of a « 0.11 cms“^, which is consis tent with accelerations ob¬ 
served in other BAL quasars (e.g.,rVilkoviskii & Irwinll200U 
[Rupke et^l2002tlGibson et aLlUOPSik 

In the standard model for BAL quasars the wind proper¬ 
ties should depend on the X-ray/UV SED and the proper¬ 
ties of the shielding gas. Indeed, C IV BALs are observed 
to be stronger and to have higher velocities in quasars with 
greater levels of X-ray absorption indicating that the shielding 
gas plays an import ant role in determining the properties of 


the outflowing wind ( 

Brandt et al.ll2000l iLaor & Brandtll2002l 

iGallagher et al.ll2006 

ilGibson et al. 1120091 iWu et alJl20i nh. In 


PG 1004H-130 the amount of X-ray absorption does not vary 
significantly with time, but despite this, large differences in 
the BAL EW are created. This further suggests that the X-ray 
absorption we observe is not due to the typical shielding gas. 
The X-ray absorption could be due to the BAL outflow itself. 
The presence of the C IV ion in the UV material requires 
an io nization parameter in the range log -0.3 to -I- 1.4 
(e.g.. iKallman & McCravl[T9^ . This is consistent with the 
values observed when the X-ray absorption is modeled with a 
partially ionized model (see Table|2]l. 

We are unable to test for any relation between the X-ray 
flux and the EW of the UV BAL as only 2 UV and X-ray 
observations were made at similar epochs. However, despite 
large differences in EW, the X-ray flux is similar in each case, 
suggesting no relationship. This is perhaps not unexpected 
as the observed X-ray emission in PG 1004H-130 is likely 
mostly from a small-scale jet (either because the continuum 
source is X-ray weak, or the X-ray normal continuum source 
is absorbed by Compton-thick shielding gas). This flux may 
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not be expected to influence the BAL production strongly. 

5.2. Nature of the Extra-nuclear X-ray Emission 

X-ray emission from extended jets has now been detected 
in over 100 AGNO but the physical processes involved in 
its production are still often uncertain. PG 1004H-130 shows 
such X-ray emission, a knot, in the same direction as the SE 
FRI radio jet; i.e., at a PA of 130° located at a distance of 8" 
(30 kpc) from the core. No X-rays are detected in relation 
to the NW FRII radio counterjet, but this is not unusual; the 
majority of detected X-ray jets are one-sided, with the X-rays 
corresponding to the brighter of the radio jets (lWorrallll2009h . 
The HYMOR radio nature is often interpreted as a difference 
in the environment on either side of the quasar; i.e., the jet 
propagating into a denser medium to the SE is quickly decolli- 
mated giving the FRI appearance, whereas the jet propagating 
to the NW through a more tenuous m edium appears as an FRII 
(e.g.. iGopal-Krishna & Wiit^ 120001) . Shocks, created when 
the je t collides with an obstacle (e.g., IB landlord & Koenigll 
Il979t) . accelerate electrons which then cool by emitting syn¬ 
chrotron X-rays and appear as bright X-ray knots. This sce¬ 
nario natur ally explains w hy knots are commonly obseryed 
in FRI jets (I Worralll 120091) as the probability of such shocks 
occurring is greater due to the higher density of the medium. 

Some X-ray knots haye shown exceptional yariability, such 
as the knot in th e jet of M87 which increased in intensity by 
50% oyer 5 yrs (iHarris et al.ll2006h . Howeyer, no eyidence 
for yariability was o bseryed in the X-ray knots in Cen A 
(iGoodger et al.|[2M^ . or in the quasar 3C 273 (l.lester et al.l 
120061) . Although the knot in PG 1004H-130 does not show sub¬ 
stantial yariability between the 2 Chandra obseryations, the 
1 .6a decrease in flux (a factor of 2 in ~ 7.5 rest-frame years) 
is consistent with electrons that haye been excited in a shock 
and are cooling yia synchrotron emission, with typical loss- 
times for such electrons being ~ tens of yrs. Expected power- 
law spectral indices for this synchrotron emission are F ^ 1.6 
(lAchterberg et al.lbOOlt) . which steepen to F 2.1 when en¬ 
ergy losses become more important. Although we do not see 
a significant difference in the power-law slope between ob¬ 
seryations, largely owing to the low number of counts in the 
spectra and hence large errors on the parameter, the yalues we 
obserye are broadly consistent with this scenario. 

Howeyer, the X-ray knot is spatially extended by 4" 
(15 kpc). Since the loss-times of the electrons are too short 
for them to diffuse across this distance before emitting X-rays, 
the emission is likely being produced from a combination of 
multiple shock sites. This may explain why the flux decrease 
(and F increase) is not highly significant, as we are obserying 
the combined properties of these sites. 

We note that in alternat iye models for the X-ray emission 
such as beamed IC/CMB (iTayecchio et al.ll200C)l: ICelotti et al.l 
120011 ) . no flux drop is expected as the typical loss-times for the 
lower-energy electrons inyolyed are ~ 10® yrs. This scenario 
is also disfayored as the jet in PG 1004H-130 is likely inclined 
at 0 > 45° (M06) implying little line-of-sight beaming. 

The peak of the X-ray knot emission lies closer to the 
core than the peak of the radio jet emission. Offsets 
such as this are also commonly o bseryed in other sources 
(e.g., iHarris & Krawczynskill200^ . and are often explained 
as a single population of electrons cooling by emitting ra¬ 
diation at progressiyely lower frequencies from X-ray to ra- 
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dip as they propagate down stream along the jet (iBai & Feel 
l2003HHardcastle et al.ll200^ . In such a scenario optical emis¬ 
sion may be expected between the X-ray and radio peaks, 
but this is not obse r yed in 7/57" imag ing of PG 1004H-130 
(iBahcall et al.lll99'^ iMiller et al.l 120061) . This perhaps indi¬ 
cates that the X-ray and radio emission instead originate from 
different electron populations. 

The deeper combined exposure time of both Chandra obser¬ 
yations allows us to further constrain the nature of the diffuse 
X-ray emission which was preyiously obseryed in M06. Fig- 
ure|5](bottom right) shows some of this emission extending to 
a radial distance of 65" (250 kpc) from the core to the SE, un¬ 
derlying the same area in which the broad FRI radio emission 
is obseryed, which includes both the jet and the lobe. Simi¬ 
lar emission is obseryed to the NW, although it only extends 
to 55" (210 kpc) and hence does not trace the location of the 
FRII lobe which is located further from the quasar. This sug¬ 
gests that the diffuse X-ray emission in PG 1004H-130 is more 
closely related to the radio jet emission, rather than the lobes. 
The possible detection of a radio counterjet on small scales 
(C. C. Cheung, 2014, priyate communication) also agrees with 
this interpretation. 

X-ray spectra of the diffuse emission, extracted from the 
green wedges shown on Figure |5] (bottom right), are well-fit 
with a power law. This suggests that non-thermal processes 
are responsible for the diffuse emission, a likely origin being 
unbeamed IC/CMB. The diffuse X-ray emission is stronger 
to the SE (Lo. 5 - 4 keV = 1-81!^;^ x lO^^ergs”') than the NW 
(Lo. 5 - 4 key = 8.56;!;3'2 X 10^'ergs”'). This is consistent with 
the distribution of colder gas implied by the HYMOR nature 
of PG 10044-130, which is likely responsible for the stronger 
radio jet to the SE and the weaker counterjet to the NW. 
Electrons produced in the IC/CMB model are of a relatiyely 
low energy, and hence haye lifetimes of ~ 10® yrs allowing 
them to trayel large distances before emitting X-rays. The 
diffuse emission could therefore haye an origin related to the 
jet and counterjet, but be located surrounding the radio jet as 
obseryed. 

6. SUMMARY 

We haye presented new Chandra, XMM-Newton, and HST 
COS obseryations of the RF BAF quasar PG 10044-130. We 
haye compared these to archiyal X-ray and UV data in order 
to inyestigate potential changes in the X-ray and UV absorp¬ 
tion, and in the X-ray jet. We also combine both Chandra 
obseryations creating a deeper exposure with which to inyes¬ 
tigate the nature of the X-ray jet and diffuse X-ray emission 
in more detail than preyiously possible. We summarize our 
main results below. 

1. The amount of intrinsic X-ray absorption present in 
PG 10044-130 does not appear to yary significantly be¬ 
tween each of the 4 Chandra and XMM-Newton ob¬ 
seryations. In each case, the absorption can be well 
modeled with a partially coyering absorber, or a par¬ 
tially ionized absorber (best-fitting spectral parameters 
are listed in Table |2|i. The low column densities. Ah = 
8 X 10^°-4 X 10^' cm“^ for a fully coyering absorber, 
and lack of yariability suggest that the absorption may 
not be due to the typical shi eldin g g as in yoked in quasar- 
wind models. See Sections [3.2l and l5n 

2. The X-ray flux of PG 10044-130 yaries by up to ~ 40% 
oyer 612 days (obseryed frame; see Figure [T]i, but re- 
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mains at least a factor of 15x weaker than is expected 
for a typical RL quasar with a similar UV and radio 
luminosity. As the source does not show significant 
variability in the amount of X-ray absorption present, 
these changes are more likely related to intrinsic varia¬ 
tions within the contin uum sou rce i tself, or fluctuations 
in the jet. See Sections (3. 1 1 and lSTI 

3. A strong C IV BAL is observed in our new 2014 
HST COS observation, with EW = 11.24 ±0.56 A. 
This is significantly variable, with an absolute change 
of IS.EW = 8.26 ±0.66 A and a fractional change of 
AfiVT /{EW) = 1.16±0.11 from the previous 2003 HST 
STIS observation, 3183 days earlier in the rest-frame of 
PG 1004H-130. While the absolute change in EW is not 
exceptional, the fractional change is on e of the high est 
observed in a BAL quasar. See Sections [3.3l and l5m 

4. An X-ray knot in the FRI radio jet is clearly detected 
in the combined, ^ 100 ks Chandra data (see Figure |5] 
top panels). It appears downstream from the peak of 
the radio emission, and at a distance of ^ 8 " (30 kpc) 
from the central X-ray source. It has a spatial extent 
of ^ 4" (15 kpc). Its spatial and spectral properties are 
consistent with synchrotron emission from populations 
of shock-heated electrons. No similar X-ray cou nter¬ 
part to t he FRII counterjet is detected. See Sections [4. II 
and l5.2l 

5. Asymmetric diffuse X-ray emission is detected around 
PG 1004H-130, with a lower luminosity underlying the 
unseen counterjet to the NW (Lo. 5 - 4 keV = 8.56^3 j x 
10"^*ergs“') than the X-ray detected jet to the SE 
(L(). 5 _ 4 kev = l-81^g44 X 10"^^ erg s“'). In both cases ather- 
mal model for the emission is disfavored and a non- 
thermal origin, possibly due to inverse Compton scat¬ 
tering of CMB photons, is more likely. The SE emis¬ 
sion has a best-fitting power-law slope of T = 1.65^° 74 , 
consistent with the value for the NW emission, T = 
1.95^} 05 . In both cases, the intensity of the emission 
decreases relatively smoothly, radially from the central 
source (see Figure | 6 ] bottom). It extends out to 65" 
(250 kpc) to the SE and 55" (210 kpc) to the NW. See 
Sections I4~2l and l5.2l 

In this work we have conducted a comprehensive, high S/N, 
multi-epoch. X-ray and UV spectral monitoring campaign of 
the RL BAL quasar PG 1004H-130. The results determined 
here serve as a template for other RL BAL quasars. The 
logical next step would be to study systematically the prop¬ 
erties of a larger sa mple of RL BALs (thus extending both 
this work and that of iMiller et al.ll20()^ iSaez et aDl201^ and 
iWelling et al.ll201^ . However, the majority of such objects 
are not economical Chandra targets due to their high z and 
low X-ray fluxes. For example, to obtain 2 spectra with ^ 500 
counts for each of the 12 snapshot targets from iMiller et al.l 
(120091) would require > 2 Ms of Chandra time. The forth¬ 
coming Af/zena observatory, with its larger effective area, will 
provide a better facility for this type of study. 
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